NTUTH-01-505B 
September 2001 



Perturbation Calculation of the Axial Anomaly 
of a Ginsparg- Wilson lattice Dirac Operator 



Ting-Wai Chiu and Tung-Han Hsieh 

Department of Physics, National Taiwan University 
Taipei, Taiwan 106, Taiwan. 
E-mail : twchiu@phys.ntu.edu.tw 



Abstract 

A recent proposal suggests that even if a Ginsparg- Wilson lattice 
Dirac operator does not possess any topological zero modes in topologically- 
nontrivial gauge backgrounds, it can reproduce correct axial anomaly for 
sufficiently smooth gauge configurations, provided that it is exponentially- 
local, doublers-free, and has correct continuum behavior. In this paper, 
we calculate the axial anomaly of this lattice Dirac operator in weak 
coupling perturbation theory, and show that it recovers the topological 
charge density in the continuum limit. 



PACS numbers: 11.15.Ha, 11.30.Rd, 11.30.Fs 



1 Introduction 



Recently, one of us ( TWC ) has constructed |] a Ginsparg- Wilson Dirac oper- 
ator which is 75-hermitian, exponentially-local, doublers-free, and has correct 
continuum behavior, but it does not possess any topological zero modes for 
topologically-nontrivial background gauge fields. This suggests that one might 
have the option to turn off the topological zero modes of a Ginsparg- Wilson 
lattice Dirac operator, without affecting its physical behaviors ( axial anomaly, 
fermion propagator, etc. ) at least for the topologically-trivial gauge sector. 
Therefore it is interesting to verify explicitly that it indeed reproduces the 
continuum axial anomaly for smooth gauge backgrounds, in the framework of 
weak coupling perturbation theory which is amenable to analytic calculations. 
The Ginsparg- Wilson lattice Dirac operator proposed in Ref. [|TJ is 

D = a~ 1 D c (l + rDc)^ 1 , r>0, (1) 

with 

D c = £7^, T» = fi?f, (2) 

f = Lsjv y /2 . ^=-e^- (3) 

\Vt 2 + w 2 + wJ ^ 

Here 7^ is the naive lattice fermion operator and — w is the Wilson term 
with a negative mass — c ( < c < 2 ) 

n^v) = \\U^x)5 x+ ^ y - Ul(y)5 x ^ y } , (4) 
Uft{x) = exp 



iagA. 



a 

X + -fi 



(5) 



r 



Gf L al + (J v (j\ = 25^ , 



(6) 
(7) 



w(x, y) = c--^2 [2S x>y - Uf 1 (x)5 x+ ^ y - Ul(y)5 x ^ )y \ , < c < 2 , (8) 

where the Dirac, color and flavor indices have been suppressed. Note that in 
(|TJ), we do not fix r = l/2c, but keep it as a free parameter, since we intend 
to show that the axial anomaly is independent of r in the continuum limit. 
Evidently, the lattice Dirac operator ([!]) is 75-hermitian 

D ] = 75^75 , (9) 
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and satisfies the Ginsparg- Wilson relation |2j 



£75 + 75^ = 2raDj 5 D . (10) 

In the free fermion limit, D is exponentially-local, doublers-free, and has 
correct continuum behavior ||. These properties should be sufficient to guar- 
antee that D can reproduce correct axial anomaly in a smooth gauge back- 
ground. 

In the next section, we calculate the axial anomaly of ([[[) in weak coupling 
perturbation theory, and show that, in the continuum limit, it recovers the 
topological charge density of the gauge background 



32tt 2 , 
where 

F„ u = d„A v - d v A„ + ig[A^ A v \ . 

Before we proceed to the perturbation calculations in the next section, we 
first clarify that even though a Ginsparg- Wilson lattice Dirac operator can 
reproduce the continuum axial anomaly for sufficiently smooth gauge config- 
urations, it does not necessarily imply that it could possess topological zero 
modes for topologically-nontrivial gauge backgrounds. This can be seen as 
follows. 

Consider a gauge configuration with positive definite topological charge 
Q > 0, and its topological charge density p(x) is very smooth for all x, 



d A x p(x) = Q . (11) 

Then, in principle, p(x) can be decomposed into two parts, 

p(x) = p w (x) + p (x) , (12) 

where p w (x) is the "winding number density" which is positive for all x, and 
it is the smoothest particular solution of (|ll|), 



d A x p w (x) = Q , Pw(x) > Vx , (13) 

and po(^) is the solution of the homogeneous equation of (|TT1), 

J d 4 x p (x) = . (14) 

If one assumes that the amount of the topological charge Q is bounded by 
a constant times the square root of the space-time volume, then in the infinite 
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volume and continuum limit, p w (x) becomes infinitesimally small at any x, 
thus it becomes unobservable. In other words, even if one calculates the axial 
anomaly A(x) (|TjJ) of a Ginsparg- Wilson lattice Dirac operator, and shows 
that it agrees with p(x) ~ Po{x) in the infinite volume and continuum limit, 
however, one still cannot be sure whether the integral / d 4 xA(x) is equal to 
the topological charge Q or not. 

In order to have a definite answer to this question, one needs to perform a 
nonperturbative calculation of the axial anomaly on a finite lattice. Otherwise, 
it could hardly detect the "winding number density" p w (x) at a site. However, 
to our knowledge, there is no preceding nonperturbative analytic calculations of 
the axial anomaly for any Ginsparg- Wilson lattice Dirac operator, on a finite 
lattice. Therefore, in this paper, we do not intend to tackle the difficult prob- 
lem of performing a nonperturbative analytic calculation of the axial anomaly 
on a finite lattice. 

Nevertheless, from the argument presented in Ref. 0, it is clear that ([!]) 
does not possess any topological zero modes. Thus, it follows that, on a finite 
lattice, the axial anomaly for at least some of the lattice sites could not agree 
with the topological charge density, since its sum over all sites is equal to 
the zero index, rather than the nonzero topological charge Q of the gauge 
background. Therefore, for any topologically-nontrivial gauge background, we 
have no difficulties to understand why the axial anomaly of ([I]) can not be 
exactly equal to the topological charge density for all sites on a finite lattice. 
However, in the infinite volume and continuum limit, p w {x) — > 0, then the axial 
anomaly A(x) can appear to agree with the topological charge density p(x) 
for sufficiently smooth gauge backgrounds, even though the integral / d 4 xA(x) 
is always zero. Thus, for a sufficiently smooth gauge background ( except for 
any possible "winding" at the infinity ), we can calculate the axial anomaly 
( at any point except the infinity ) of the lattice Dirac operator (|I]) in weak 
coupling perturbation theory, and show that in the continuum limit, it agrees 
with the topological charge density of the gauge background, regardless of the 
global topological charge. 



2 Perturbation Calculation of the Axial Anomaly 

The axial anomalyQ of a flavor-singlet of Ginsparg- Wilson lattice fermions can 
be written as j3| 

A(x) = tr[75(I — raD) (x, x)] 

(2^/(2^ el{P ~ q) ' X ^75(1-^)0^)] , (15) 



1 It is understood that the axial anomaly (15) at a lattice site has to be normalized by 
the volume of the unit cell, a 4 , when it is compared with the topological charge density in 
the continuum. 
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where the trace runs over the Dirac, color and flavor space. Now we expand 
D in power series of the gauge coupling ( through the link variables ) 

D = A) + D 1 + D 2 + D 3 + D 4 + 0(g 5 ) (16) 

where D n denotes the terms containing the factor g n , the gauge coupling to 
the n-th power. Then we have 

I-raD 

= e - e rfie - e (d 2 - d^d^eo - e (d 3 - d 1 e d 2 - d 2 e d 1 + die Q die di)eQ 
— eo(c?4 — d\eod 3 — d^e^di — d 2 eod 2 + dieodieod 2 + dieod 2 eod\ + 

+d 2 e die di - dieod^odieod^eo + 0(g b ) (17) 

where 

e = l-raD = (1 + rfo)" 1 , (18) 

di = r£7,(/itf/o + /otf/o + /o«{>/i), (19) 
n 

d 2 = r J2 iMh + /otf /o + /otf/a + fittfo + /itf/i + /otf/i) , (20) 
n 

ds = ^7,EE(/4/3-H)» (21) 
d4 = r^l^iZUi^h-i-j) , (22) 

At i=0 3=0 

/ = /0 + /l + / 2 + /3 + /4 + O((? 5 ) , (23) 

t" = t£ + t!t + 1% + 1% + 1% + 0{g 5 ) . (24) 

Evidently some terms in ([H]) do not contribute to the axial anomaly, since 
they only have a factor of a product less than four distinct 7 matrices, thus 
they vanish after taking the trace with 75 in (|l^). Note that each term of d n 
only has one factor of 7 matrix, and eo contributes at most one factor of 7 
matrix [ see Eq. ( [33|) ]. Thus, the terms e , and e d n e ( for any n ) can be 
dropped from the series, when ( |T?D is substituted into (|T5|). 

Then the axial anomaly at the order g 2 is due to the fourth term in (|17"D, 

A(aO = / / / e^>MlMp)di(p,k)e (k)d 1 (k,q)e (q)] (25) 

Jk 

where 

d 4 p 



(27T) 



4 " 



etc. 



Similarly, substituting higher order terms of (|17]) into dl5[) gives the ax- 
ial anomaly at orders g 3 and g 4 , respectively. Although there seems to be 
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many terms involving higher order vertices ( d 2 ,ds, etc. ) at higher or- 
ders, the nonvanishing contributions at orders g 3 and g A only come from 
terms e§d\e§d\e§d\e§ and eodie die dieodie respectively, as we shall show 



later. At this point, it is instructive to identify eo = (1 + do)~ fllg) with 
the free fermion propagator, d\ flT9"|) with the " quark-quark-gluon" vertex, 
c?2 p0| ) with the " quark-quark-gluon-gluon" vertex, etc. Then the expres- 
sion tr(7 5 eo<iieo<iieo) contributing to the axial anomaly at 0(g 2 ), (p5|), cor- 
responds to the triangle diagram coupling the axial current to two external 
gauge bosons through an internal fermion loop. Similarly, the expressions 
t r (75eo^ieo^ieo^ieo) and tr (7560^1 eoG?ieo<iieo<iieo), contributing to the axial 
anomaly at 0(g 3 ) and 0(g 4 ), correspond to the quadrilateral ( box ) and 
the pentagon diagrams coupling the axial current to three and four external 
gauge bosons, respectively. These are exactly the Feynman diagrams^] which 
contribute to the axial anomaly in continuum nonabelian gauge theories [H. 
Note that even though higher order vertices are allowed at finite lattice spac- 
ing, the symmetry forbids their contributions to the axial anomaly, thus only 
the "quark-quark-gluon" vertex d\ can enter the fermion loop contributing to 
the axial anomaly. 

In the following three subsections, we evaluate the axial anomaly at orders 
g 2 , g 3 and g A respectively. 

2.1 The Axial Anomaly at the Order g 2 

To evaluate (pip, first, we obtain d\ flT5| ) by expanding / and in power series 



of the gauge coupling, (p3|) and fl24|) 
Consider 



h = Vt 2 + w 2 + w = h + h + h 2 + h 3 + h 4 + 0(g 5 ) , (26) 
which satisfies the identity / • / • h = 2c, i.e., 

C/o + fi + ■ ■ -)(/o + fi + ---)(h + h 1 + ---) = 2c. (27) 
To order g, it gives 

fofoho = 2c 
fifoh + fofiho + fofoh = . 

Solving these two equations, we obtain 

1/2 

fo( P , k) = I r= 2 4= ) — k) = f (p)6 4 (p - k), (28) 



Itlip) + wl{p) + w (p) 

- /o 2 (p)/o 2 (fc) 
2cf (p) + f (kY 



fiip,k) = ~ /^ ) f S ^ 1 M {p,k), (29) 



2 See, for example, Figure 22.2 in Ref. 
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where 

tl( P ) = -£tf(P)#(p) = £sin 2 (p M a) . (30) 
Similarly, we obtain 

h , 7 \ ™l(P, fc)Mg) + Wp(fc)] ~ tf(p, fc)[tg(p) + tff(fc)] 

^(p)+^o 2 (p) + ^( fc )+ w o 2 W 

(31) 

Here we recall some well-known basic formulas in weak coupling perturba- 
tion theory ( see Appendix A ), 

P(p,k) = t^p)5 A (p-k)+t1(p,k)+0(g 2 ) , 
w(p,k) = w (p)5 4 (p-k)+w 1 (p,k)+O(g 2 ) , 



where 



*o(p) = «sin(p M a) , 
Wo(p) = c - £[1 - cos(^a)] 



%(p,k) = gA„(p-k)dX 
wi(p,k) = ^gAftip-fydnWa 



' p + k" 



2 y 
a 



Mp-k) = J2^ i(p ~ kHx+mA ^[ x + ^ 

Therefore, the axial anomaly (|25|) at the order g 2 can be written as 
A 2 (x) = g 2 (I e^" 1 ^ Y.I H^ip-k)A v {k-q)}x 

J v J a ,, J k 



where 

eo(p) = 1 - raD (p) 



tT[^ 5 eo(p)di,n(p,k)eo(k)di jU (k,q)e (q)] , (32) 
1 1 -rf$(p)-yt (p) 



l + d (p) l + r 2 /o 4 (^o(p) 

= 6(p) + £ 7 ^(p), (33) 

do(p) = r/ 2 (p) 7 -t (p)=E7^(p) , (34) 

di(p,fc) = gJ2A»(p-k)d ltlx (p,k) , (35) 



di^ip, k) = r[/i j(U (p, fe)7 ■ to (AO /o (AO + /o(p)7 ■ *o(p)/i,^(p, fc ) 

+ /o(p)7 M ^( £ ^)/o(fc)l , (36) 



/lrflCP ' } " 2c /o (p) + /o(fc) X 

n /p + ifc\ (2?) Mp) + w (k)) - dX (*±*) (tf (p) + (*)) 



v / ^(p) + ^o 2 b) + v /t oW + ^) 



(37) 



In the limit k = p, (|37|) and ( p6|) reduce to 



/i>,p) = fy/o(p), (38) 
diAViV) = ^[ r /o(p)7"to(p)/o(p)] = <VoG») • (39) 



To evaluate the integral in fl32|), we change the variables p — > p + k and 
g — > g + k, then ([32]) becomes 



A 2 {x) =g 2 e i( v-«>* Y,^{Mp)M~Q)} G»„(p,q) (40) 



pjq 



where 

G^v{p,q)= j tT[j 5 e Q (p + k)di^(p + k,k)eo(k)di jU (k,q + k)e (q + k)] . 
In the following, we show that 

d d 

G^{p,q) = a a G ^ fa g) lp,g=o + °( a ) 

A (T CPA Cg^ 

= -M(c)^4p A g ffVAff + 0( a ) (41) 

which, when substituted into (fffiP, leads to the axial anomaly at the order g 2 , 
A 2 (x) = g 2 M(c) J2 ^xMi^A v - d v A»){d x A a - d a A x )(x)} , (42) 

where M(c) is a coefficient which tends to ^2 ( for < c < 2 ) in the 
continuum limit. 

We expand G^p, q) in power series of p and q, 

d d 
G^fag) = G^{0,0) + Y,P\-^—G tJ , u (p,q)\ m=0 + Y, ( l<T-^—G fl v(p,q)\p, q= o 

X °P^ a (7(7 

d d 

+ ^2PxQa ~ a ^Wfa g) Ip,9=0 + °( Q ) • ( 43 ) 
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It is easy to see that the zeroth order and the first order terms in (^) vanish, 
by symmetry and the basic properties of 7 matrices, 

tr(Ts) = tr( 757 Q ) = tr( T67 V) = tr^VY) = tr(7 5 7V7V7 <5 ) = 0, (44) 
tr( 75 7 a 7 /3 7 A 7 CT ) = tr(l)e Q/3A(T . (45) 

Explicitly, 



G^(0,0) = / tr{ 75 e (k) d^d (k) e (k) d u d (k) e (k)} 

Jk 

= tr(l) / e ap5a c a (k)c\k) b(k) dX(k) d u d°(k) (46) 



where (|33|), Q3§, (|3|), (@|) and (gD have been used. Note that repeated 
indices are summed over in Eqs. fl46|) , (|50"D , (|5l"D, ([54]) and (|55|). Evidently the 
integrand in (|46| ) vanishes due to contraction of the completely antisymmetric 
tensor e a p$ a with the symmetric tensor c a (k)c 5 (k), i.e., 

Y,^(3S,c a (k)c\k) = 0. (47) 

aS 

Hence, we have 

G^(0,0) = . 

Next consider the first order term of the power series (pl|), 
d 

tr[ 75 e (/c) dpd 1 ^(p,k)\ p=k e (k) d ljV (k,k) e {k)) 
+ tr[ 75 d x e (k) d 1>t ,(k, k) e (k) d h „(k, k) e (k)\ } (48) 
From (p6|), it is easy to see that 

d^(p,k)\ p=k = Y,rhl^ (49) 

where each term has one factor of 7 matrix, and the explicit expression of 
is not required for our purpose. 

Then the first term and the second term of the integrand ( f48|) both vanish, 

tr[ 75 e {k) 7 • h Xfl (k) e (k) d v d Q (k) e (k)} 
= tr(l) e Q/M<7 c Q (A;)c 5 (A;) 6(fe) h^(k) d v d°{k) = , (50) 



and 



tr[ 75 d x e (k) d^d Q (k) e (k) d u d (k) e ] 

tr(l) (e^ + e^) cP{k)b{k) d x c 5 {k) d„d%(k) d u d° (k) 

+ tr(l) e a/3(j5 ^(k)c s (k) d A &(A;) S^W 9 Mk) 

, (51) 
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where (0), (0), (f§), (g|), (g|) and (0) have been used. 
Thus 

GV(p,g)| Pi9 =o = . 



Similarly, we have 



<} G M1/ (p,g)U = o = . (52) 



Now consider the second order term of the power series fl43|). 

n — t; — Gfj, u (p, q)\ pq= Q 

tr[7 5 d x e {k) d hlx (k,k) e {k) dgd hv (k,q)\ q=k e (k)] 

+ tr[7 5 e (k) <9*d liM (p, k)\ p=k e (k) d 1;V (k,k) d a e (k)} 

+ tr[7 5 e (k) dpd Ul (p, k)\ p=k e (k) d^d 1 ^(k,q)\ q=k e (k)] 

+ tr[7 5 d x e (k) di tlx (k,k) e (k) d ltU (k,k) d a e (k)]} . (53) 

Evidently, the first three terms of the integrand are zero, by symmetry. Ex- 
plicitly, the first term is equal to 

tr[7 S d x e (k) d^d (k) e (k) 7 • h av (k) e (k)] 
= tr(I) (e 5a ^ + e 5ai(3 ) c?(k) b(k) d x c s (k) d^(k) h^k) 

+ tr(!) e Q ^ 5 c?{k) c\k) d x b(k) d^(k) hl v {k) 
= 0, (54) 

the second term is also zero since it has the same form of the first, and the 
third term is 

tr[ 75 e (k) 7 ■ h Xfl (k) e {k) 7 ■ K v (k) e {k)) 
= tr(I) e aPSl c a (k)c s (k) b(k) h^(k) hUk) = , (55) 

where (0), (0), ([§), (||), (@), (§|) and (0) have been used. 
Thus, (|53|) becomes 

9 -G^(p,g)| M=0 



tr{9 A [7 5 e (fc)] d^do(k) e (k) d v d (k) <9 CT e (£;)} . (56) 



From ( j33|) and ([34]), it is easy to see that e (k) also satisfies the Ginsparg- 
Wilson relation 

75 + l {k)^e {k) = 2^e Q {k) , (57) 



9 



which, after differentiation with respect to k\, gives 

[dxeo \k)]>y 5 eo(k) + e ~ 1 (k)j 5 [d x e (k)} = 2d x [ l5 e (k)} • (58) 

Substituting (|58|) and cfo(&) = e^" (A;) — f [ from Eq. (|i~8l) ] into the integrand 
of fl56|), then the integrand becomes 

-tr{d x eo 1 (k) 75 e (A;) d^ l {k) e (k) d v e^ l {k) d a e (k)} 
+ ^tT{e \k) 75 <9 A e (A;) d^(k) e (k) d u eQ l {k) d a e Q (k)} 
which can be further reduced to 

-- tr{7 5 d^e (k) d u e^ l (k) d a e (k) dxe^ik)} 
+\ tr{ 75 d x e (k) d^\k) d u e (k) d a e^(k)} , (59) 
where the identities 

e (k) d fl eo 1 (k) e (k) = -d^e (k) (60) 
e \k) d^k) eo\k) = -d^\k) (61) 

have been used. By symmetry, it is obvious that the contribution of the second 
expression of (^) is equal to that of the first one, thus fl56|) becomes 

IfiiyXa = Tj q G ^ {j>, q) \ p q=0 

d 4 k 

tr{7 5 d^e Q {k) d v e^{k) d x e (k) d a e^(k)} 

= /l^^rr^) ^^tt^) a * rf ° w } (62) 

where the domain of integration is the 4-torus, T<± = ®| =1 [— 7r/a, vr/a], in which 
the endpoints ( ±7r/a ) in each direction are identified to be the same point. 
In the limit a — > 0, T 4 — > ®f =1 (— oo, oo), which is invariant under the trans- 
formation 

^^-^T> // = l,---,4, (63) 

for any r^O. 

Now we change the variables according to (0), then (|62|) becomes 



/ (H tr (rrhm) %m fcdbo ) (64) 
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where the arguments of do are = — {r 2 k^)~ l . Note that dk^d^ is invariant 
under the transformation (|63D. 

Therefore, our ansatz of evaluating ( 161 ) in the continuum limit is to sub- 
stitute d (K) by dQ l {k) in the integrand of (|64D, since d (k) — > ira^ as a — > 0, 
and the trace with 75 picks out four distinct 7 matrices, one from each of the 
four factors with partial derivatives. Then (|64|) reads as 



tr 1 i , . 7T\ (fc)9 A : I , 7T\ (fc) 



d 4 A; 



rtr 



{75 cU" 1 ^) <9 A /i(*0 ^^(fc)} , (65) 



where 



(2tt 

d (k) 



h(k) 



d (k) 



The integral (|65"D can be evaluated by first removing an infinitesimal ball 
B e of radius e from the origin ( k = ) of the 4-torus Zi, then performing the 
integration, and finally taking e — > 0, i.e., 

W = T7T7 lim / rf4fc tr i75 ^A(fc) d./T 1 ^) d A /i(fc) ^/T 1 ^)} (66) 
= T^lim/ d 4 fc 9 M tr{ 75 d v hT\k) d x h(k) d a hr\k)} (67) 

167T 4 e^O Jr 4 \B e 

where the operation in ( CT ) produces (BH), plus three terms which are 
symmetric in /iz/, fi\, and fia, respectively, hence neither one of these three 
terms contributes to the integral. 

Then according to the Gauss theorem, the volume integral over T^\B e can 
be expressed as a surface integral on the surface S t of the ball B t , provided 
that the integrand is continuous in T4 \ B e . The last condition is satisfied 
since h(k) = raD (k) is analytic, and h{k) 7^ ( free of species doublings, for 
< c < 2 ) for any k e T 4 \ B t (|. Thus (BTJ) becomes 



W = lim / d3s n Ml5 h(k) d v h~\k) d x h(k) d a h-\k)} (68) 

167T 4 e-»0 J5 e 

where n M is the /i-th component of the outward normal vector on the surface 
S e . Since g?o(&) - ► irafi as — > 0, we can set do(^) = irafc on the surface 
and obtain 

t 1 r f j3 A ( h \ lv ( 7a \ 7/ 



16tt 4 ^o75 £ " l' d \l + r 2 a 2 A; 2 M 2 Vl+r 2 a 2 Py P 
16vr 4 " v " y ^ >A,J ^0 75^ " £; 4 (l + r 2 a 2 P) 



' tr(I) lim / d 3 s — (69) 
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where we have used (fH|), (0), and the property 

f d 3 s n il k v F{k 2 ) = 5^ j d 3 s n il k il F{k 2 
Finally, the result of the integral (^9|) is 



IS 








/ dc 


b / d6 2 sin 6 2 


10 


Jo Jo 



I --^Le Mm 1 

\xv\a - 47r 4 C ^A CT am (l + r2a2e2)2 

where iV/ denotes the number of fermion flavors. Note that (|70|) is independent 
of the parameter r in D ([!]). From ( [7D|) and (41), we obtain 

M ( c ) = ^' 0<c<2, (71) 
and the axial anomaly (|42| ) at the order g 2 , 



^ 2 (x) = f-^ E e ^ tr{(9^ - d v A„){d x A a - d a A x )(x)} , (72) 



where 

= / e tp x A,(p) , etc. 
J P 

2.2 The Axial Anomaly at the Order g 3 

Inserting the order g 3 terms of (|17|) into (|15|), we getf] 

A 3 (x)=g 3 f [e i ^- x {J a (p,q) + J b (p,q)}-g 3 [ f [ e l ^ *J c (p,r,q) (73) 

Jp Jq Jp Jr J q 

where 

<^ab,g) = / tr{7 5 e (p)rfi(p, k)e (k)d 2 (k,q)e {q)} , 

■^(p,?) = tT{'y 5 eo{p)d 2 {p,k)eo{k)d l (k,q)eo(q)} , 

Jk 

J c (p,r,q) = tT{~f 5 e (p)d 1 (p,r)e (r)d 1 (r,k)e (k)d 1 (k,q)e (q)} , 



and eo and c?i are given in (|33|) and (|35|) respectively. It is easy to show that J a 
and Jb vanish in the continuum limit, and only J c has nonzero contribution to 
the axial anomaly at the order g 3 . To proceed, we derive d 2 (|20| ) by expanding 
/ and in power series of the gauge coupling, (1231) and (E^), similar to the 



Hhe term e^d^e® has been dropped since tr^eoc^eo) = 0. 
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procedure of deriving d\ in Section 2.1. First, we solve for fa from the higher 
order equations in (|27"D , then together with the formulas of t£ ( |126| ) and w n 
( |127| ) derived in the Appendix A, we obtain 

^(p, ^) = 5 2 ^ if (p - fc)d 2 ,(p, fc) = 9 2 ^ - ^(p, fc) , (74) 

/Lt fall 

where Ay{jp — k) is defined in ( |118| ), and the explicit expression of d^Jp, k) 
is not required for our subsequent calculations. In general, to any order of g, 
we have 

i(p I *) = s n E^(p-*)^(p. fc ) = 9 n E^ ) (p-*KC(p> fc ) > ( 75 ) 

fl faU 

which can be used for higher order calculations. 

Using (|35]) and (|74]) , and changing variables p —>■ p + k and q — > q + k in 
J a and J^, while p ^ p + r + k, r — > r + k and g — > q + in J c , then we can 
rewrite ([F^) as 

A 3 (x)=g 3 [ fe^-^{I a (p,q) + I b (p,q)}-g :i [ [ [ e^ +r -^I c (p,r,q) (76) 

Jp J q Jp J r J q 



where 



I a (p,q) = £tr[i>)i( 2 )(-g)]Gr(p,<7) , (77) 
/ 6 (p,g) = Y,^\ P )A u (-q)}Gr(p,q) , (78) 
/ c (p,r,g) = E tr K(p)^W^A(-9)]Gr A (p,r,g) , (79) 

and 

GT(p, q) = I tr[7 5 e (p + k)d lifl (p + fc, k)e (k)d 2 , v (k, q + fc)e (g + k)] (80) 

Gb U {p,q) = / tr[7 5 e (p + k)d 2 ^{p + k, k)e (k)d 1)V (k,q+ k)e {q + k)} (81) 

G^ uX (p, r,q) = [ tr[7 5 e (p + r + k)d ltfl (p + r + k,r + k)e (r + k) x 

d ltV (r + k, k)e {k)d ltX {k, q + k)e {q + k)} . (82) 



k 



In the following, we show that G a and vanish in the continuum limit. 
First, we expand (^0) in powers series of p and g, 

GT(P,9) = G^(0,0) + ^p^G^ a %p,0)\ p=0 + ^:q^G^(0,q)\ q=Q 
+0(a) , (83) 
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where the zeroth order and first order terms can be easily shown to be zero, 
by symmetry. Explicitly, 



GT(0,0) = / trfr 5 e (k)diAk, k)e (k)d 2 ,„(k, k)e (k)\ = , 

Jk 

since the integrand vanishes identically, 

tr[7 5 e (/c)9 / ,e " 1 (A;)eo(A;)(i 2)! ,(fc, k)e (k)) = -ti[-f b d fl e {k)d 2 ^{k,k)e {k)] = 



where Eqs. (g3j), (^9]), O and (g§ have been used. 

To evaluate the first order terms in (^), we observe thatQ 



d 



d 



J2l X h X ^(k) = 7 " K,{k) 



dq 



-d 2 ,u(k,q + k)\ q=0 = J2^9au( k ) = 1 ■ 9au{k) 



(84) 
(85) 



A 



which follow from fllTf), fl20|), ( |36"D and flT4|), where the explicit expressions of 



h x ^{k) and g^ v {k) are not required for our subsequent calculations. Then 




GfM)| 



p=0 



tr{75[<9<Te (/i:)<9^ (£;) + e (fc) 7 • h a ^(k)}e (k)d 2 ^(k, k)e (k)} = , 

since its integrand is equal to 

tr(l)e^ A(5 {[9 (J 6(A ; )^(A;) + 6(A;)^(A;)]c /3 (A;)c 5 (A;)4 l/ (A;)} 
+ tv(l)(e apxs + e am ){[d a c a (k)dX(k) + c a {k)h^(k)]c\k)d 5 2 ^k)} = , 
where (^), (@), (H), (@)> (fH) and (flTD have been used. Similarly, we have 

1^(0, g)U 

tr[^eo(k)d^d {k)eQ{k)Y) ■ g au (k)e (k) + d 2:U (k, k)d a e (k)}} = . 

Therefore G% u (p,q) = 0(a) — > in the limit a — ► 0. By the same token, 
G^ u (p, q) = in the continuum limit. 

Next we expand G c in power series of p, q and r, 

GfV r, g) = GT A (0, 0, 0) + £ j^G^fa 0, 0)| p=0 

+ E^^r A (0,r,0)| r=0 +X:? CT ^^(0,0,g)| 9=0 + O(a) . (86) 
4 Note that each term of di (|l9|) and c?2 (|20|) has only one factor of 7 matrix. 
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It is easy to see that the zeroth order term vanishes, 
Gf A (0,0,0) 

tr[jseo(k)d ljti (k, k)e {k)d hv (k, k)e (k)d ljX (k, k)e (k)} = , (87) 

k 

since its integrand vanishes, 

tr[7 5 eo(A;)di, M (A;, k)e Q {k)d liV {k, k)e (k)d ljX (k, k)e {k)} 
= tr[y 5 eo(k)d fl eo 1 (k)e (k)d 1/ eo l (k)eo(k)d x eo 1 (k)e (k)] 
= tr[7 5 a / ,eo(/c)^eo 1 (/c)aAe (/c)] = , 



where (|39|), (|60|) and (Q) have been used. 
Next we evaluate the first order terms 

= -^tr{7 5 [9 <T e (/i;)9 At e 1 (A;) + e (fc)7 • h afi (k)]e (k)d u eo l (k)d x e (k)] 
= ~ J k ti{ 1 sd u e {k)d^\k)e {k)d^\k)d x e {k)} , (88) 

where (|39"D, (]84]) and (|60|) have been used, and the last equality is due to the 
vanishing of the second term ( containing the factor 7-/i - /Lt (&) ) in the integrand, 

tr{7 5 e (A;)7 • h afl (k)e (k)d u eQ 1 (k)d x eo(k)} 
= ti(l)e a ^ s d x b(k)c a (k)c^h^(k)d u d 5 (k) 

+ tr(l)(e^ 7<5 + e 0a ^)b(k)c a (k)d x c s (k)h^(k)d v dl(k) 
= . 



Comparing (p8|) with (|56|), one immediately obtains^ 

^Gf A (p, 0, 0)| p=0 = l, vXa = -^^a. , (89) 



which has been evaluated in (|70|). 
Similarly, we obtain 



5g -GT A (0, 0, g)| g=0 = = ^e^ A(T . (90) 



'Note that 9 M e 1 (fc) = d^doik), from (33) 
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Next we evaluate 
d 



Gr A (0,r,0)| ? 



^{tr[75<9 CT eo(/c)9 M e " 1 (A;)eo(A;)<9 ! ,eo 1 (A;)eo(A;)9Ae " 1 (A;)eo(A;)] 

+ti[-f 5 e (k)d^e 1 (k)d a eo(k)d u e 1 (k)e (k)d x e 1 (k)e (k)] 
+tr[<y 5 eo(k)d a d l _ l d (k)e (k)d 1/ eQ 1 (k)e (k)d x e 1 (k)e (k)} 
+ tr[ 75 e (A;)^eo 1 (A;)e (A;)(7-^(A;))eo(A ; )9 A e 1 (A;) eo (A;)]} (91) 
, (92) 



where in the integrand of (pTf) , the first two terms cancel each other, and the 
last two terms vanish respectively, after the ( by now familiar ) manipulations 
using fl), (H), (0), @, (0), © and ©. 



Substituting fl87|), (H), © and (H) into (H), we obtain 

GfV r, q) = --4 J>, - q a ) e, uX(T + 0(a) . (93) 
Therefore, in the continuum limit, the axial anomaly (|76|) at 0(g 3 ) is 
M%) = HTirE^/ I [ e^ +r -«>%p a - q„MA„(p)A v (r)A x (-q)) 

^ [xv\a J P Jr J 1 

9 3 N f 



■ V) -i E e^xatr^idaA^x) - d^A a (x))[A u (x), A x (x)} + 



[A„(x),A„(x)](d a A x (x) - d x A a (x))} , (94) 



where 



A,{x) = / e ip - x A^(p) , etc. 



2.3 The Axial Anomaly up to Order g 4 

Inserting the order g A terms of (|TT| ) into Ql5"D, we havef] 

A(x) = 9 4 JJ^-^{j 1 (p,q) + J 2 (p,q) + J 3 (p,q)} 

-9 4 JJJ q e l{p - qyx {Mp, r, q) + J 5 (p, r, g) + J 6 (p, r, g)} 
+9 4 I I I I e^-^J 7 (p } r,s,q) (95) 

Jp Jr Js Jq 

6 the term eo^eo has been dropped since tr^eocLieo) = 0. 
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where 

Ji{p,q) = / te[j5eo(p)di(p,k)e (k)d 3 (k,q)e (q)} 
Jk 

J2(p,q) = / tr[7 5 e (p)cf 2 (p, k)e {k)d 2 {k,q)e (q)} 

Jk 

h{p,q) = / tr[7 5 e (j9)rf 3 (p, k)e (k)d 1 (k,q)e (q)} 
Jk 

J^{p,r,q) = / tr[7 5 e (p)di(p,r)e (r)di(r, k)e {k)d 2 {k,q)e {q)} 

Jk 

Jh(jP,r,q) = / tr[7 5 eo(p)di(p,r)eo(r)d 2 (r, k)e (k)d l (k,q)e (q)] 

Jk 

Je(p,r,q) = / tr[7 5 e (p)d 2 (p, r)e (r)4(r, k)e (k)di(k, q)e {q)} 

Jk 

J 7 (p,r,s,q) = / tT[j 5 e (p)d 1 (p,r)e Q (r)d 1 (r,s)e (s)d 1 (s,k)e (k)d 1 (k,q)e (q)] 

Jk 

Now we change the variables as follows : 

(i) p — > p + k and q — > q + k in Ji, J 2 , and J 3 ; 

(ii) p^p + r + /c, r — > r + fc, and q — > g + k in J 4 , J 5 , and J 6 ; 

(hi) p-^p + r + s + A;, r— >r + s + fc, s — > s + fc, and q — > g + in J 7 . 

Then fl9~5|) becomes 
AW = 9 4 JJ q e^-^{l 1 (p,q) + I 2 (p,q) + h(p,q)} 

-9 4 JJJ q e l(p+r - q) x {h(p, r, q) + I 5 (p, r, g) + J 6 (p, r, q)} 
+9 4 I I I I e l ^ r+s -^I 7 (p,r,s,q) (96) 

Jr Js 

where 

h(p,q) = Y.^\Mp)Af\-q)]GT{p,q) 

h(p,q) = T,^\p)Af\- q )}G^{p iq ) 

hip A) = T, t ^ ) (p)A u (~q)]G^(p,q) 

hfar,q) = Y, t 4Mp)Mr)A?(-q)}Gr X (p,r J q) 

h(p,r,q) = Y,^[Mp)^ ) (r)A x (-q)]Gr\p,r,q) 

h(p,r,q) = Y, t AA? ) (p)Mr)A x (-q)]G^\p,r,q) 

I 7 (p,r,s,q) = £ tx[A^p)A v {r)A x {s)A a (-q)]G^{p,r,s,q) 
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and 

Gi U (p, q)= tr[7 5 e (p + k)d hfM (p + k, k)e (k)d 3tU (k, q + k)e (q + k)] 

J k 

G 2 u (p, q) = tr[7 5 e (p + k)d 2 ^(p + k, k)e (k)d 2 , v (k, q + k)e (q + k)] 

J k 

G 3 v (p, q) = / tr[7 5 e (p + k)d 3 ^(p + k, k)e (k)di tV (k, q + k)e (q + k)] 

J k 

GT X (p, r,q) = / tr[j 5 e (p + r + k)di^(p + r + k,r + k)e (r + k) x 

J k 

d l)V (r + k, k)e (k)d 2t \{k, q + k)e (q + k)] 
G^ uX (p, r,q) = / tr[j 5 e (p + r + k)d\^{p + r + k,r + k)e (r + k) x 

Jk 

d 2 , u {r + k, k)e (k)d hX (k, q + k)e (q + k)} 
Gq vX {p, r,q) — / tr[7 5 e (p + r + k)d 2 ^(p + r + k,r + k)e (r + k) x 

Jk 

di,„(r + k, k)eo(k)di t \(k, q + k)e (q + k)] 
G? uXa (p, r, s,q) = / tr[7 5 e (p + r + s + k)d 1<fl (p + r + s + k,r + s + k) x 

J k 

e (r + s + k)di jV (r + s + k, s + k)e (s + k) x 
di jA (s + k, k)e (k)d 1>a (k, q + k)e (q + k)} 

where (p5|), (|74l) and (|75|) ( for n = 3 ) have been used, and iW(p), etc. are 
defined in Eq. (|TT8| ). 

In the following, we show that all G's vanish in the continuum limit except 
G 7 , by expanding each one of them in power series of p, q, r and s, respectively. 

First consider 

GT(P, q) = GT(0, 0) + 0{a) , (97) 

where 



GT(0,0) = / tr[7 5 eo(fc)di, M (fc,fc)eo(fc)d v (fc,fc)eo(fc)] = , 

Jk 

since its integrand vanishes, 

tr[jse (k)d lyfl (k, k)e (k)d 3tU (k, k)e {k)) 
= ti{l)e aM& b(k)c a (k)c\k)dXmiAk) = , 

where (|3"3"D, fl39|), ([T5|), (f44|), (|45| ) and p7| ) have been used. Similarly, we obtain 

Gf (0,0) = Gf (0,0) =0 . 

Next consider 

Gf A (p, r, g) = Gf A (0, 0, 0) + O(a) , (98) 
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where 

Gf A (0,0,0) 

tv[j 5 e (k)di !tl (k, k)e (k)d 1>v (k, k)e (k)d 2 ,x(k, k)e (k)} = , (99) 

since its integrand vanishes, 

tr[y 5 eo( y k)d^eo 1 (k)eo(k)d I/ do(k)eo(k)d2,x(k, k)e {k)} 
= -tT[-f 5 d fl eo(k)d u d ( y k)eo(k)d 2t x(k, k)e (k)} 
= -tr(I) {e a ^ 5 d,b(ky(k)c\k)d u d«(k)dl x (k)} + 

+tr(I){(e a/37 i + e a/357 )6(A ; )^(fc)^c Q (fc)^^(A;)4 iA (A ; )} 

= , (100) 

where ([^), ©, © and (0) have been used. By the same 

token, we have 

Gf A (0,0,0) = GT A (0,0,0) =0 . (101) 
Finally, we consider 

Gf A >, r, s, q) = Gf A<T (0, 0, 0, 0) + 0(a) (102) 

where 

Gf Act (0,0,0,0) 

tr[7 5 e (£;)di^(fc, k)e (k)d ljV (k, k)e {k)d ljX {k, k)e (k)d 1>t7 (k, k)e {k)} 
tr[7 5 e (A;)9 A1 eo 1 (A;)eo(fc)<9 I ,eo 1 (/i;)eo(/i;)9 A e 1 (A;)eo(/i;)9 (T eo 1 (A;)eo(/i;)] 
ti[ lb dM k ) d u%\k)e (k)d x e Q \k)d a e^)} • (103) 



Here we have used (|39| ) and (pOf) to simplify the integrand. Comparing ( |103|) 
with (|56|), one immediately obtains 

Gf A,7 (0, 0, 0, 0) = I^ Xa = -^e^ uXa , (104) 

which has been evaluated in (|70|) . 

Thus, in the continuum limit, only G7 has nonzero contribution to the 
axial anomaly. Substituting (|102|) into (p6|) , we obtain the axial anomaly at 
the order g A , 

M?) = -jr E e ^ CT / / / / ^^(pAwWh 

4 



p J r J s J q 



TV 

-^-^ XI e^A CT tr{yl M (x)^(x)A A (x)A (T (x)} 



8tt 



o 4 ^ 



327T 2 

(JLUA(T 



]T e^Aatr{[A /l (x), 4,(x)] [A A (x), A CT (x)]} . (105) 
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Finally, adding (|72|), ( p^ ) and ( |105| ) together gives the axial anomaly ( |i5|) 
in the continuum limit, 



A(x) = tr [75 (1 — arD)(x, x)] 
a 2 N f 

= %d E ^ tr{i^(x)F A(7 (x)} , (106) 

where 

i^, = 9^ - d v A^ + A„] . (107) 

Note that the axial anomaly (|106|) is invariant for any r > 0, since the anomaly 
coefficient ( JT0| ) is independent of r. 

This completes our perturbation calculation of the axial anomaly of the 
Ginsparg- Wilson lattice Dirac operator proposed in Ref. 



3 Concluding Remarks 

Several remarks are as follows. 

It is obvious that (|106|) also holds for other (H) such that D is doublers- 



free, exponentially-local, and has correct continuum behavior. In particular, 
(|106T) holds for 

f=(^TT=¥—-) > °< c < 2 > ( 108 ) 

\ \Jt 2 + W 2 + Wj 2 

as proposed in Ref. IJ. 

It is instructive to examine how well the continuum axial anomaly can be 
recovered on a finite lattice, by evaluating the integral /1234 (0) as a numerical 
sum over the discrete momenta on a finite lattice. In Table 1, the ratios of 
A234 to 1/87T 2 are listed for several lattice sizes ( L 4 ), as well as for a range 
of r, respectively. Here the lattice spacing a is set to one, and the value of 
c [ see Eqs. (Q) and (g) ] is fixed at 1.0. Evidently, the integral /1234 tends 
to the continuum value 1/87T 2 as L — > 00, independent of the parameter r. 
This also provides a verification of our ansatz to substitute d (K) by dQ l {k) 
in evaluating I^ v \ a ( |64| ) in the continuum limit. 

Next we repeat the same calculations of /1234 for another /, namely, re- 
placing (0) by 

f = /, 2 + 2 V • (109) 

The results are listed in Table 2, which show that the integral /1234 ( as a 
function of lattice size L 4 ) approaches the continuum value 1/87T 2 much faster 
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L 


r 




0.5 


1.0 


2.0 


4.0 


16 


0.9204 


0.9800 


0.9848 


0.7617 


32 


0.9810 


0.9952 


0.9988 


0.9907 


64 


0.9953 


0.9988 


0.9997 


0.9999 


128 


0.9988 


0.9997 


0.9999 


1.0000 


256 


0.9997 


0.9999 


1.0000 


1.0000 



Table 1: The ratio of the integral I1234 [ Eq. (|62 ) ] to I/Stt , for lattices of 
sizes 16 4 ,32 4 ,64 4 , 128 4 ,256 4 , as well as for r = 0.5,1.0,2.0,4.0, respectively. 
The lattice spacing a is set to one, and the parameter c [ Eqs. (§) and (||) ] is 
fixed at 1.0. 



L 


r 




0.5 


1.0 


2.0 


4.0 


16 


1.0000 


1.0000 


0.9911 


0.7565 


32 


1.0000 


1.0000 


1.0000 


0.9914 


64 


1.0000 


1.0000 


1.0000 


1.0000 


128 


1.0000 


1.0000 


1.0000 


1.0000 


256 


1.0000 


1.0000 


1.0000 


1.0000 



Table 2: The ratio of the integral I 12U [ Eq. (§2j) ] to 1/8tt 2 , for / = 

2c/ (v^ 2 + w 2 + w). Other parameters are the same as those in Table 1. 

than that in Table 1. This indicates that (|109| ) may be a better choice than 
(H), especially for small lattices. 

In general, (0) may not be in the form Tjj = ft^f, then our present 
perturbation calculation may not go through without modifications. We refer 
to our former derivation || for the general case, though some of our inter- 
mediate steps need further clarifications. If one assumes that the coefficient 
of the axial anomaly, g 2 /327r 2 , is the same for U(l) and SU(n) background 
gauge fields, then it can also be determined || by imposing a gauge configura- 
tion with constant field tensors to Luscher's formula @ for the axial anomaly 
of Ginsparg- Wilson lattice Dirac fermions in a U(l) background gauge field, 
provided that D is topologically-proper. However, for a topologically-trivial 
Ginsparg- Wilson lattice Dirac operator such as (|1|), it seems to be necessary 
to perform an explicit calculation in order to determine its axial anomaly. In 
passing, we also refer to other axial anomaly calculations |], [II], [XT], for 
the overlap Dirac operator |13|, |T4]| , as well as that in the original Ginsparg- 
Wilson paper M. So far, a nonperturbative analytic calculation of the axial 
anomaly for any Ginsparg- Wilson lattice Dirac operator, on a finite lattice, is 
still lacking. 
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In summary, we have shown that the lattice Dirac operator (|1|) reproduces 
continuum axial anomaly for smooth gauge configurations, even though it does 
not possess any topological zero modes in topologically-nontrivial gauge back- 
grounds. If one insists that the topologically zero modes of a lattice Dirac 
operator are crucial for lattice QCD to reproduce the low energy hadron phe- 
nomenology, then one should assure that a Ginsparg- Wilson lattice Dirac op- 
erator is indeed topologically -proper, before it could be employed for lattice 
QCD computations. However, so far, there does not seem to have compelling 
experimental evidence that these topological zero modes are physically rele- 
vant, unlike the axial anomaly in the trivial gauge sector, which accounts for 
the decay rate of the neutral pion. So there might be a very slight possibility 
that lattice QCD with topologically-trivial quarks could reproduce low energy 
hadron phenomenology. These issues seem to deserve further studies. 

A Basic Formulas in Weak Coupling Pertur- 
bation Theory 

In this appendix, we set up our notations by deriving some basic formulas in 
weak coupling perturbation theory. These formulas are used in our anomaly 
calculation in Section 2. 

Consider the forward and backward difference operators 




(110) 




where the link variables ( for SU (n) gauge group ) are defined as 





(113) 



Expanding V 1 * 1 



in power series of the gauge coupling g, 




(114) 



and performing the Fourier transform 



^e-^-^)(v;) n (x,y) = (v;) n ( P ,g) , 



(115) 



we obtain 



(Vj) (p, q) = (e*" a - 1) 5\p -q) = (V+)„(p) 5\p - q) 



(116) 
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and 



x,;y 



n! 



-i(p-g)-(a;+oA/2) _ 



5 



x+afi,y 



TV. 



: {iag) n 



A* ^ M ' 



where 



(V+) 



A^Qr) — g) — e -i(p-g)-(g+oA/2) 
P + *A _i(p M + 9/i )a/2 _ j _ 



A„ ( X + -fb 



W n M = ^(p-^(v;) (^) 



(117) 

(118) 
(119) 



Similarly, we have 

(V;) (p,g) = (l-e-*".) < 5*(p_ g ) = (v;) o (p)5*(p- g ), (120) 



(121) 



Then we obtain the weak coupling perturbation series of t M (f|) and w 
in the momentum space : 



where 



w(p, q) 



w (p) 

KM 

w n (p,q) 



toW\p -q) + tf(p, q) + #(p, q) + .. 

c "^E[V+(p,g) + V M (p,g)] 
z n 

Wo{p)5 i (p -q) + wt{p, q) + w 2 (p, q) + 



isin(p M a) , 

c - ^[1 - cos(p M a)] , 



^E4 n) (p-#> 



(122) 
(123) 

(124) 
(125) 

(126) 
(127) 
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